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Abstract: Multiple inflation is a model based on N=l supergravity wherein there are 
sudden changes in the mass of the inflaton because it couples to 'flat direction' scalar fields 
which undergo symmetry breaking phase transitions as the universe cools. The resulting 
brief violations of slow-roll evolution generate a non-gaussian signal which we find to be 
oscillatory and yielding /nl ~ 5 — 20. This is potentially detectable by e.g. Planck but 
would require new bispectrum estimators to do so. We also derive a model-independent 
result relating the period of oscillations of a phase transition during inflation to the period 
of oscillations in the primordial curvature perturbations generated by the inflaton. 



Keywords: Cosmic Microwave Background Radiation: non-gaussianity, Early 


Universe: inflation, supersymmetry and cosmology 





1. Introduction 



A major goal in modern observational cosmology is to detect any non-gaussianity of the 
temperature anisotropies in the cosmic microwave background (CMB) M, B H. All single 
field, slow roll, inflation models give rise to anisotropies that are gaussian Therefore 
detection of a significant non-gaussian signal would falsify such models and provide new 
insights into the dynamics of inflation. 

Quantifying such a signal is however not a straightforward task. Whilst gaussianity is a 
well-defined property, non-gaussianity is not and the anisotropies can, in principle, deviate 
from gaussianity in many different ways |J. Also, the CMB temperature anisotropies 
measured by WMAP are very close to gaussian ||. Therefore any measure designed to 
detect non-gaussianity needs to be sensitive to a very small signal. 

One measure of non-gaussianity that is particularly useful is the three point correlation 
function, or 'bispectrum' of the temperature anisotropies. Gaussian random variables (and 
functionals of gaussian random variables) have the property that all odd power correlation 
functions are zero. This makes the bispectrum the lowest order statistic for which any non- 
zero result would indicate a departure from gaussianity. The bispectrum contains much 
more information than the power spectrum as, in general, it depends on both scale and 
shape. Therefore, if a non-zero bispectrum is detected it will also be an extremely useful 
statistic for constraining models of the early universe. 

However there are limitations to how much information can be inferred from the bis- 
pectrum: 

• Modern CMB experiments possess very large numbers of pixels of data. e.g. N ~ 10 6 
for WMAP and ./V ~ 10 7 for Planck ||. Considering that each higher order of correla- 
tion function will require an additional factor of N calculations, exact determination 
of the bispectrum quickly becomes impossible. Therefore various estimators have 
been constructed 0, but each can look only for a specific type of bispectrum and so 
might miss a signal different from the type being searched for. 

• Higher order correlation functions suffer more from cosmic variance because more 
information is required at each scale to compute them. In this work we will deal with 
a scale dependent signal, therefore limitations due to cosmic variance are particularly 
relevant. 

These problems are not necessarily insurmountable. Computing power will only con- 
tinue to get stronger, making bispectrum estimators increasingly powerful with time. For 
each potential signal, estimators can also be constructed for that specific signal, ensuring 
processing time is used as efficiently as possible. The estimator can be optimised for de- 
tecting primordial non-gaussianity while discriminating against secondary non-gaussianities 
arising from e.g. unsubtracted point sources or residuals from component separation ||. 
Finally, in models that have non-trivial scale dependence, a correlation will likely exist 
between the power spectrum and bispectrum as we illustrate in this paper. 
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One of the first methods for quantifying non-gaussianity involved rewriting the New- 
tonian potential as ]|], 

= $ L (x) + /niMW - (^(x)}), (1-1) 

where /nl is a constant, and <&l(x) is a gaussian variable. However, this parameterisation 
has a rather specific form and does not capture other possible deviations from gaussianity. 
A more general parameterisation can be obtained by allowing /nl to depend explicitly on 
the wavevector, k. In terms of the adiabatic curvature perturbation, the parameterisa- 
tion becomes || , 

C(x) = Cl(x) - ^/nl * (Cf» - (Cf»»- (1-2) 

Here, the ★ product is used because /nl has been allowed to depend on scale. 1 The factor 
of -3/5 comes from the relationship between the adiabatic curvature and the Newtonian 
potential at matter domination. 

For /nl as defined in Eq. (|i~l~l) , the WMAP 5- year constraint is -9 < /nl < HI H- 
This reinforces the point made already: the primordial temperature anisotropics are very 
close to gaussian. Given that the amplitude of these anisotropies is AT/T ~ 10 -5 the 
non-gaussian contribution to these anisotropies is at most ~ 10~ 8 (or ~ 0.1% of the overall 
anisotropy). In general the expectation from inflation is that the anisotropies should be 
close to gaussian, therefore this result can be considered as a success of the inflationary 
paradigm. 

More precisely the expectation from inflation of non-gaussianity as parameterised by 
Eq.(|1.2j) depends on the specific model. It is known that for single field, slow-roll inflation 
with a canonical kinetic term and a vacuum initial state, the result is /nl ~ e *C 1 j§] where 
e is the usual slow-roll parameter defined later in Eq. (|3.6| ). The best sensitivity expected 
from the Planck satellite is to /nl of 0(5) while the secondary contribution from post- 



inflationary evolution is of 0(1) [10]. Therefore the prediction of the simplest inflationary 
toy models 2 is that there should not be a detection of primordial non-gaussianity in the 
near future. 

Conversely a deviation from the simplest toy models can produce a larger value for 
/nl- The study of these effects serves two purposes. Firstly, from the cosmological perspec- 
tive, a detection of /nl will immediately rule out single field slow-roll inflation and focus 
attention on determining how inflation actually occured. Secondly, from the perspective 
of inflationary model building based on fundamental physics, as the bounds on /nl grow 
tighter, some interesting models can be ruled out if there is no detection. 

Much work has been done on non-gaussianity generated due to multiple scalar fields 
(e.g. Ref.Jll]]), non-canonical kinetic terms (e.g. Ref . |i~2|] ) or non- vacuum initial states 
(e.g. Ref. fill ) . However, surprisingly little attention has been paid to the possibility of 
non-gaussianity generated by a violation of slow-roll. Ref. [14] did consider this in a context 



1 See Eq.(229) in Ref. M for how to define this product. 

2 By this we mean a generic fine-tuned potential such as the frequently used V((f>) = m 2 (f> 2 which has 
not yet been convincingly obtained from a physical theory, especially since inflation occurs in such models 
a t 4>> M P . 
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when multiple fields are present; however, the non-gaussianity itself is generated by the 
multiple fields, not the violation of slow-roll (see also Ref . |l5| ) , Refs.|l6|, 17] do consider 



the non-gaussianity generated by violation of slow-roll, but for a toy inflationary model. 
We follow their methods closely but consider a physical model of inflation. In both models 
the violation of slow-roll generates sharp features in the power spectrum and also generates 
a ringing in the bispectrum with a characteristic period which we calculate below. 

1.1 Features in the Power Spectrum of Primordial Fluctuations 

It is an expectation in the simplest toy models of inflation that the power spectrum of 



temperature anisotropies in the CMB should be nearly scale-invariant [18]. However this 
need not be the case for physical models. If the observed spectrum was indeed perfectly 
scale-invariant, such models would be ruled out and naively this might seem to be the 
case. Assuming the "concordance" ACDM cosmology, the primordial spectrum is well 
parameterised as a power law with spectral index n s = 0.960 ± 0.014 ||. This indicates 
that there cannot be any significant scale dependence that affects the spectrum over a 
large range of scales but it does not preclude the existence of e.g. localised oscillations 
in the spectrum, or other sharp features. Moreover since the observed anistrotopies arise 
from the convolution of the unknown primordial spectrum with the transfer function of 
the assumed cosmological model whose parameters are being determined, it is obvious that 
both unknowns cannot be determined simultaneously without further assumptions. 

In fact, there are indications that the primordial spectrum might not be a scale-free 



power law, even assuming the ACDM cosmology |19|, |2Cj, |21|, |22j, [23|, |24J[. At large angular 
scales (> 50°) there is essentially no power and there are anomalous 'glitches', especially in 
the range of multipoles I ~ 20 — 40 [^, pq| . The statistical significance of these anomalies 
is not sufficient to claim a definite detection, however it is strong enough to provide some 
tension with the fit to a scale-free power-law primordial spectrum which has only a 3% 



probability of being a good description of the WMAP 1-year data [25], although this did 
improve to ~ 7% with the WMAP 3-year data release |26|]. Future measurements of the 
EE and TE mode polarisations by the Planck and (proposed) CMBPol satellites will throw 



light on whether these glitches are real or not [27|. 

"In the absence of an established theoretical framework in which to interpret these 
glitches (beyond the Gaussian, random phase paradigm), they will likely remain curiosi- 
ties" [ p6[ . Indeed if there were a model that had predicted, without ambiguity, the po- 
sition and amplitude of the glitches, this would be seen as very strong evidence for the 
model. Although no such model exists, the general possibility of generating glitches over 
a range of scales had in fact been proposed prior to the WMAP observations in the con- 
text of 'multiple inflation' wherein the mass of the inflaton field undergoes sudden changes 
during inflation f28| . This generates characteristic localized oscillations in the spectrum, 
as was demonstrated numerically in a toy model of a inflationary potential with a 'kink' 
parameterised as f29fl: 



V{4>) = \m 2 c\> 2 



1 + c tanh 



(1.3) 
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By tuning the position, amplitude and gradient of the kink, the locations of the glitches 
can be varied to match the glitches seen in the power spectrum. One can then perform 
statistical likelihood tests to determine whether the fit is better with the glitches, but 
with the additional parameters, or with the simple scale-free spectrum. It was found by 
the WMAP team that the fit to the 1-year data improves significantly (by A% 2 = 10) 
for the model parameters <j) s = 15.5 Mp, c = 9.1 x 10~ 4 and d = 1.4 x 10~ 2 Mp, where 
M P = (8ttGn) _1/2 ^ 2.44 x 10 18 GeV pi. This analysis was repeated later using the 



WMAP 3-year data, with similar results [ 31 1 . 

This seems encouraging, however m in the toy model above is not the mass of the 
inflaton — in fact in all such monomial 'chaotic' inflation models with V oc (j) n , inflation 
occurs at field values cj)\ n ^ > Mp, hence the leading term in a Taylor expansion of the 
potential around <j)- m ^ is always linear in (J) (since this is not a point of symmetry), rather 
than quadratic as for a mass term [32]. The effect of a change in the inflaton mass can be 
sensibly modelled only in 'new' inflation where inflation occurs at field values c/>i n fl << Mp\ 
and an effective field theory description of the inflaton potential is possible. The 'slow-roll' 
conditions are violated when the inflaton mass changes due to its (gravitational) coupling to 
'flat direction' fields which undergo thermal phase transitions as the universe cools during 
inflation |psfl . The resulting effect on the spectrum of the curvature perturbation was found 
by analytic solution of the governing equations to correspond to a 'step' followed by rapidly 
damped oscillations [j33|. 3 

The next step should be to predict other observables, having used the power spectrum 
to constrain all the parameters in the model. In this paper we calculate the bispectrum 
of the multiple inflation model [28], using its predicted power spectrum [33] and the set 
of parameters which provide the best reduced \ 2 m the ACDM cosmology [35]. We also 
examine the effect on the bispectrum of varying these parameters over their full natural 
range. 4 The non-gaussianity is found to be potentially detectable by the Planck satellite, 
or perhaps even a reanalysis of the WMAP data. 



2. Multiple inflation 

The biggest difficulty in inflationary model building is obtaining a potential that is both 
sufficiently flat to support inflation for ~ 60 e-folds and stable towards radiative corrections. 
N = 1 supergravity (SUGRA), the locally realised version of supersymmetry (SUSY), is a 
natural framework for achieving this (see Ref.[]36| for a comprehensive review, especially of 
the cosmological issues discussed below). In SUSY there are usually many 'flat directions' 
in field space, i.e. scalar degrees of freedom which, while SUSY remains unbroken, have 
perfectly flat potentials. When SUSY is broken, the flat directions are lifted and usually 
acquire a mass-squared related to the SUSY breaking scale. It would seem natural to 
identify one of these flat directions as the inflaton but to achieve a sufficient number of e- 
folds of inflation, the mass-squared of the inflaton needs to be much smaller than the Hubble 

3 Although a similar phenomenon had been noted earlier for the case where the inflaton potential has a 
jump in its slope |$4j, such a discontinuity has no physical interpretation. 

4 We use the word 'natural' in this context to mean "stable towards radiative corrections" . 



- 4 - 



expansion rate during inflation — if it is not, the evolution is too quick and inflation ends 
after only a few e-folds. However, due to SUSY breaking by the vacuum energy driving 
inflation, the natural expectation is mh ~ H 2 . This is commonly known as the r\ problem 
in supergravity and superstring model building because it relates to the slow-roll parameter 
r\ (defined in Eq.|3.7|) being too large to support inflation. 

There has been much work done to attempt to circumvent the rj problem (see Ref. 



for a review and Ref. [38, |39j for some recent attempts in the framework of string/brane 
models). In this work we will assume that the inflaton mass has a symmetry protecting it 
against SUSY-breaking corrections. This can be done for example by having the inflaton 



in a 'hidden sector' where it interacts with the other sectors only gravitationally [40|. 

In the simplest inflationary models, scalar fields other than the inflaton itself are usually 
ignored. This is justified by the argument that, because of their larger masses, they will 
not contribute to the adiabatic density perturbations produced during inflation. If these 
fields decay before the end of inflation then their decay products will be diluted by the 
expansion of the universe and they will also not contribute any isocurvature perturbations. 
Such fields usually cannot couple to the inflaton either because any such coupling could 
endanger the flatness of the inflaton potential. 



Multiple inflation [28] is based on N = 1 SUGRA and assumes that while 'flat direction' 
scalar fields have their natural SUSY-breaking scale masses of order the Hubble parameter, 
there is just one scalar field (the inflaton) with a mass that is kept anomalously small 
through a symmetry. The amplitude of the potential during inflation will be the SUSY 
breaking scale (which need not be the electroweak scale as SUSY breaking can be different 
during and after inflation). The flat directions couple gravitationally to the inflaton and 
thus affect its evolution when their own masses change as the universe cools during inflation 
and they undergo symmetry-breaking phase transitions. 

2.1 Evolution of the flat direction fields 



It is assumed [28] that before inflation begins the universe is in a thermal state at tem- 
perature T. At this point the flat directions, tp, are trapped at the origin by a thermal 
potential oc tp 2 T 2 . Unless protected by a symmetry, the mass \i of the field will be of the 
order of the Hubble parameter, H, which is itself determined by the scale Vo ~ A 4 of the 
vacuum energy driving inflation. Henceforth we work in units where the reduced Planck 
mass Mp = 1, therefore fi 2 ~ H 2 and H 2 ~ A 2 /3. 

Although the mass-squared is negative, the flat directions will be lifted at large field 
values by non-renormalisable operators oc ip n (Mp = 1). For small field values, at non-zero 
temperatures, the full potential for the flat directions is thus |28| ], 

V(tp) = V + (~ Y + uT 2 ^ ijj 2 + yip n . (2.1) 

Before inflation, when the temperature is much larger than the mass of the field, it is 
trapped at the origin by a thermal barrier (quantum tunneling through which is negligible). 
When inflation starts, the temperature drops rapidly and after ~ ln(Mp/A) ~ 10 e-folds 
of inflationary expansion the thermal barrier disappears. The field can then evolve to 
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its minimum at £ = (/i 2 /n^) l l n ~ 2 and meanwhile there are ~ ln(£Mp/A 2 ) e- folds more 
of inflation [28]. 5 The field evolves (as a critically damped oscillator) according to the 
equation of motion: 

dV 



i>\ (2.2) 



when it reaches the minimum the oscillations of the field are damped rapidly (within 1— 
2 cycles) by the 3Hip term (see Fig.l of p3]). Damping due to particle creation is far 



less effective since this is happening in an inflating background [28]. Consequently no 
isocurvature pertubations are generated by such oscillations. 

The relationship of the flat directions to the infiaton 

It was shown in Ref.|28| that a term K^fitp 2 in the Kahler potential with k of 0(1) is 
consistent with the underlying symmetries and corresponds to a term ^(f) 2 ip 2 in t ne scalar 
potential where A = kH 2 . The full potential, including the infiaton, is thus: 

V(<f>, ^) = V - \m 2 ^ 2 - Y, Q/^i " \\<Ptf ~ 7Vf ) • (2-3) 

There will be a running of the infiaton mass too, however the field value of the infiaton 
remains small throughout inflation so we can neglect the higher-dimensional operators that 
would parameterise such running. 

Hence as each field, tpi, rolls to its minimum at £j, the effective infiaton mass-squared 
changes by AjS 2 . The amplitude of the potential also changes as each field falls into its 
minimum, however, this change will be very small compared to the dominant term in the 
potential V f2§. 

During periods of slow-roll, the amplitude of the primordial perturbation spectrum 
can be expressed in terms of the potential and its derivatives as, 

H 2 V 3 

Pk = s^T = ^v^' {2A) 

Where the * indicates that the quantities should be evaluated as the relevant scales 'cross 
the horizon', i.e. when k 2 = 2a 2 H 2 . A transition in the mass of the infiaton, though not 
causing a significant change in V, will cause a significant change in V' . If Aj is positive, the 
mass will decrease causing a jump in the amplitude of the spectrum, while if it is negative 
then the amplitude will fall. As the flat direction oscillates in its minimum, the mass of 
the coupled infiaton will briefly oscillate as well. If the amplitude of the oscillations is large 
enough, one would expect a 'ringing' in the primordial power spectrum [35] for which there 



is tentative observational indication |21], 22, E3, U, 



The WMAP measurements indicate that, from I ~ 40 through to the limit of I ~ 
800 set by signal-to-noise considerations, the inferred primordial spectrum is nearly flat. 
However, as was noted earlier, there is tentative evidence of a step at very large scales of 



5 The field evolves as ip oc e Ht so the occupation number of thermal states drops as ~ e~ e . There is 
no non-gaussianity generated due to the thermal fluctuations while the field is at the origin (c./. E3]) since 
observable perturbations today leave the horizon much later. 
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order the present Hubble radius [T^, 21, 22]. It was in fact a prediction of the original paper 



on multiple inflation p|| that there should be of 0(1) such feature in the ~ 10 e-folds of 
inflationary history probed by the CMB. 

The advantage of the multiple inflation model over the kink potential (EHl) considered 



in Ref.[16, [2S|] (or the toy model considered in Ref . ]34j ) is its grounding in a consistent, 
particle physics framework. While the latter model is able to improve the fit to the glitches 
in the data p0| , |3l| ] the constraints derived on the model parameters do not tell us anything 
about the relevant physical processes at these high energy scales. By contrast, constraints 
on multiple inflation are related directly to the masses and couplings of fields potentially 
present during inflation. For example, the scalar potential of the Minimal Supersymmetric 
Standard Model (MSSM) is flat along many directions in field space and a catalogue exists 
of all these directions p^j . These would be lifted during inflation due to SUSY-breaking 
by the large vacuum energy present and conceivably undergo symmetry breaking phase 
transitions. 

3. Method 

3.1 Single field justification 

We follow closely the method of Ref. [16], who use the framework of Ref.Q in calculating 
the bispectrum. This method uses the comoving gauge where, in single field inflation, the 
inflaton fluctuations are zero at all orders. When there are other, more massive, fields 
present this is no longer true because these extra fields will slightly perturb the instan- 
taneous, adiabatic, direction in field space away from the direction of the main inflaton 
field. In multiple field scenarios a simpler gauge to describe the evolution of perturbations 
within the horizon is the uniform curvature gauge. In this gauge, neglecting gravitational 
waves, the slicing of spacetime into equal time hypersurfaces is made such that the spatial 
curvature on each of these surfaces is zero. 

It was shown in Ref.Q] (and generalised beyond linear order in Ref. |p5[| ) that there 
exists a gauge- independent, non-perturbative quantity £ which, on sufficiently large scales, 
coincides with the adiabatic, scalar density perturbation at the perturbative level. Defining 
the number of e-folds of expansion as N = J Hdt, it was shown in Ref.|44|, |2!| that if we 
start in the uniform curvature gauge, this quantity can be expressed as ( = dN. This 
allows us (on sufficiently large scales) to express the adiabatic curvature perturbation as 

, ^dN^r v-^1 d 2 N .,r.,7 , olN 



2. 

I T i.j 



If we set the initial, uniform curvature, hypersurface a few e-folds of expansion after the 
modes cross the horizon then due to the much smaller mass of the inflaton, we have <p <C ip, 
hence dN/d(p (= H/cf) 3> dN/dip for all the flat direction fields. For all fields involved (in- 
cluding the inflaton), the higher order terms in Eq. fl3.1| ) are also very small. 6 Under these 



6 These terms become significant in models where non-gaussianity is seeded by entropy perturbations, 
after the perturbations have left the horizon |jq |. 
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conditions, the result £ = ^Scp holds to a very good approximation. The flat direction 
fields do not affect the adiabatic curvature perturbation and thus isocurvature perturba- 
tions do not arise gravitationally during inflation. Therefore the inflaton fluctuations will 
be zero in the comoving gauge during inflation and we are free to use the action derived in 
Ref.Q] using this gauge. 

It should be stressed however that the flat directions do play a crucial role in the 
generation of the non-gaussianities. The point of the preceding argument is that they do 
this only through their coupling to the inflaton and not directly through their coupling to 
gravity. What this means for the calculation is that when solving the zeroth order equations 
of motion for the inflaton we must include the effects of the flat direction fields. However, 
when computing the perturbative action to quadratic and cubic order in the perturbations, 
we need consider only the inflaton. 

In principle it might have been less confusing to work initially in the uniform curvature 
gauge, where the contribution from all the fields is unambiguous, and to convert to the 
comoving gauge once the curvature perturbation has left the horizon. This would have 
the benefit of less ambiguity in the present case and be a necessity in genuine multi-field 
models. In future work we intend to generalise the method of Ref.[16|, as used in this 



paper, to inflationary models involving several, dynamically equally important, fields. 
3.2 Calculating the bispectrum 

We are interested in calculating the bispectrum (three-point correlation function) of the 
quantum observable £, corresponding to the scalar, adiabatic curvature perturbation. The 
model we are considering has the standard, purely gaussian vacuum state, and the quadratic 
action of the field, £, is effectively just the free field action (i.e. a kinetic term and a mass 
term). A free gaussian field will remain gaussian, therefore in order to probe any non- 
gaussianity in the fluctuations it is necessary to consider the cubic action. We know from 
observations of the power spectrum that |C| ~ 10~ 5 i.e. higher terms in the perturbative 
expansion of the action will be less important in general. 7 Therefore, we need only consider 
the action up to cubic order. 

We proceed by re- writing the quantum operator in the interaction picture: 

G(t) =e iHt{t - to) ( 1 (t )e- iHt{t - to \ (3.2) 

where Hf is the free field Hamiltonian derivable from the quadratic action and Ci(*o) is the 
initial vacuum state value of the operator. The full quantum operator is then, 

C(t) = U 1 (t,t )\i(t)Ui(t,to), (3-3) 

where, 

U I (t,t ) = Te- lf to Hlit ' )dt ', (3.4) 

is the time evolution operator for the interaction terms (and T denotes the time ordering 
operator). 



7 This need not be true always, e.g. if a symmetry forces the cubic action to be anomalously small then 
the quartic action will produce the dominant contribution to any non-gaussianity 0]. 
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Putting these terms into the three-point function gives, 

(C(k!)C(k 2 )c(k 3 )) = (ulbitfUj) 

"'<[Ci(i) 3 ,#i(0]}dt'. (3.5) 



to 

Where the last line of equality holds to first order in the expansion of the exponentials. 

It is usually emphasised at this point that the vacuum with respect to which the 
expectation value is being taken is the fully interacting vacuum, not the free vacuum. This 
is important because the Ci appearing in these equations will be expressed in terms of the 
ladder operators that annihilate the free vacuum. Fortunately, this is not an issue in this 
particular case because, similarly to what is done in Minkowski space, one can slightly 
deform the integral contour into Euclidean space. This projects out the free vacuum, as 
explicitly shown in Ref.||. 

Up to field redefinitions proportional to the free field equations of motion, the Hamilto- 
nian in Eq.( |3.5| ) can be obtained from the action derived in Ref.[|J] (Eq.(3.9) in that paper 
and Eq.(3.4) in Ref . [|i~6f ) . This action can be expressed as a sum of terms whose coefficients 
are functions of the slow-roll parameters, defined as: 

e - ^, (3.6) 

^7h^J^ + 1P> (3 - 7) 
where the approximate equality holds when both e and r\ are small. Note that there are 
other common definitions for e and rj in terms of derivatives of the potential: ey = ^(Y' /V) 2 
and r]y = V" /V; the ones defined above are related through ey = e and rjy = —\r\ + 2e. 
These equalities hold in the slow-roll approximation when all four quantities are small. 
In terms of these quantities, the cubic action for C is:, 

S 3 = fd*x ah\C 2 + ae\{dC? - 2a\Q(dQ(d X ) + 



+^(d()(d x )d\ + a ^{dx)(dx? + 2/(C) ^ 



2 dt 

(31 



where x = e <9 C- 

The term 2/(C) § 



is proportional to the free field equations of motion 



and 



in Ref.[y is removed by a field redefinition. This can have an important effect on the 
bispectrum, however, in the present case it does not, due to its relative size. There is just 
one term in /(C) (Eq.(3.10) in Ref. [Q]) which does not include a derivative on one of the C 
terms. As we can choose to apply the correction due to the redefinition after the relevant 
scales have left the horizon, only this term will have a non-zero effect — its coefficient is 
r//4, which is much smaller than our leading term. 8 

8 Note that although r\ does become large during the phase transition in multiple inflation, we can choose 
to apply the correction after the phase transition is over, at which point 77 is again 1. The redefined 
field will evolve outside the horizon while rj is significant so this can be done only if we also calculate the 
evolution of the redefined bispectrum out to this time, which we have indeed done. 
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With the field redefinition taken care of, the Hamiltonian in Eq. (|3.q ) just becomes 
Hj = —L3. To use this to evaluate Eq.( |3.5[ ) we need to solve the commutator in this 
equation. We are free to use the standard commutation relations for (j; to do so we first 
need to express Cl(k, t) in terms of the ladder operators of the free field vacuum 

Ci(k, t) = u k (t)a k + nl k (t)a f _ k , (3.9) 

with the usual commutation relations, [a k , a^,] = (27r) 3 <5 3 (k— k') and mode functions u k (t). 
The relationship between Ci(k) and Cf( x ) * s the usual 

/d 3 k 
— C i(k)e' ikx (31Q) 



When these definitions are substituted back into Eq.(3.5), the result is proportional to 
an integral over time of the free field mode functions u k (i), evaluated on the interaction 
Hamiltonian, H\ ( [3.11| ). 

The dominant term in multiple inflation 

It has been stated already that multiple inflation violates slow-roll, nevertheless e remains 
<C 1 throughout multiple inflation for all parameter values considered here. 9 There is only 
one term in the action (3^) that does not have a factor of at least e 2 in it — this is the 
term with dr]/dt (note that x is °f order e) . However although e remains small, due to 
the transition in the inflaton mass, r/ and dr]/dt can temporarily become large. The size of 
the latter term is dictated by both the magnitude and the rate of the mass change in the 
phase transition. Multiple inflation is a small field model, hence e is very small (< 10~ 10 ) 
and the dr]/dt term dominates over all others because it has one less factor of e. 

3.2.1 Numerical method for calculating bispectrum 

Upon substituting this leading term into Eq. (|3.5| ), using the definitions of £(k) and C(x) 
and working through the commutation, the result is (in conformal time dr = dt/a): 

(C(k 1 ,r)C(k 2) r)C(k 3 ,r)) = -21m j u kl (r)u k2 (r)u k3 (r) jT dr' a 2 e|^) 3 5 3 (l>) x 

f d-y* (-[ n * \ "| ^1 







This should be compared to Eq.(3.21) in Ref.|16]. 



To obtain the bispectrum for multiple inflation, we must solve this integral. To do this 
we first need to numerically solve for the free field mode functions u k (r), the scale factor 
a(r) and the slow roll parameters e and dry/dr. Once this has been done, the results are 
put into the above integral and evaluated out to a value of r after the end of inflation. All 
that is needed after that is a rescaling of the full bispectrum to obtain a scale dependent 
generalisation of /nl- 



9 If e does become large, inflation will stop; this is potentially possible if there is a later inflationary 
epoch |28] but we do not consider this possibility here. 
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Fortunately, the slow roll parameters in the above integral are composite variables of 
other variables more fundamental to the model (e.g. e = cp 2 /2H 2 ). Therefore, we need 
only solve for these more fundamental variables and then evaluate e and drj/dt from them. 
The full list of variables that need to be solved for, with their equations of motion, are [49]: 



z" 



4 + [k ~ — K = 0, (3.12) 



with z = a\/2e, Vk = zuk, and: 

„ 2 dV 2a'<f>' 

9 +a——-\ = 0, (3.13) 

dcp a 

,.. 2 dV 2a'ib' , . 

+a 2 — + ^- = 0, 3.14 
dip a 

a" + ^ 2 + V>' 2 )-^ = 0. (3.15) 

In all of the above, the ' indicates differentiation with respect to conformal time. These 
four equations are just the free field equation of motion for the observable the two zeroth 
order equations of motion for <j) and if), and finally the Friedmann equation for the Hubble 
parameter a/a. 

We solve these equations numerically, using the Matlab function odell3 |]5C|]. For 
initial conditions we take the Bunch-Davies vacuum state jljj. This amounts to: 



Vk(To) = \hrr, (3.16) 




v' k (r ) = -i^. (3.17) 

The initial value of r = To in the integrand was set such that integration of the mode 
functions, v^, begins when the following condition is first satisfied: 

7" 

10 4 - > k 2 . 

z 

At earlier times, the mode functions are highly oscillatory and would cancel in the integral 
in Eq.( 3.1l| ). In fact there is only a small window in which this integral is non-zero — very 



early on, the mode functions are highly oscillatory and any contribution to this integral 
will be washed out, while very late on, the mode functions are frozen and thus the duk/dr 
term will be zero. Finally, when the phase transition is not occurring, the dij/dr term 
will be negligible. It is only those modes that are leaving the horizon during the phase 
transition that will give a a non-zero. contribution to the bispectrum. 

A regularisation scheme is used in Ref.[^] to counter a sharp cutoff in the integral 
in Eq.( |3.1ll ) (and improved upon in Ref . |l7[| ) . We find that a regularisation scheme is 



useful for calculational efficiency but not necessary for accuracy. It was argued in Ref. [16] 
that although the integral is technically convergent as r — > —00, this sharp cutoff would 
add a spurious contribution to the bispectrum, due to the highly oscillatory nature of 
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the mode functions at early times. We find that this is a problem only when the cutoff 
also co-incides with the violation of slow-roll. If this occurs, there was indeed a spurious 
contribution due to the sharp nature of the cutoff; however, if we begin our integration 
earlier, when the slow-roll contribution to the integral, edrj/dr is negligible, no spurious 
contribution is obtained. 

In ordinary slow-roll, to accurately calculate the contribution from the other, dominant 
terms in Eq. ([[!]) by numerical means, a regularisation scheme would be necessary. For 
multiple inflation (and for the kink model considered in Ref.|16|]) these other terms are small 
enough to be ignored. This is especially true in multiple inflation where e is exceptionally 
small. What the regularisation scheme does do is allow the integration to begin later 
without sacrificing accuracy. Nevertheless, we did not use a regularisation in calculating 
our results and altered the integration start time when there was a risk of a spurious 
contribution. 

3.2.2 Comparing a bispectrum to /nl 

We have defined a scale-dependent generalisation of the /nl parameter in Eq. fll.2p , Now we 
discuss how to calculate /nl from a given bispectrum. For a flat bispectrum a convention 
has been well established, however there is an ambiguity in the present case due to the 
scale dependence of both the power spectrum and the bispectrum. For a flat spectrum, 
following Fiefs. [||, |9|, we start by re-expressing the bispectrum as, 

(C(k 1 )C(k 2 )C(k 3 )) = (2nfd 3 P>> ) t^^. (3.18) 



In this notation, we can rewrite /nl in the form 10 



/-=-y?- <3 - i9) 

However, as pointed out in Ref.|l(|, this definition entangles some of the ringing of the 
power spectrum in the definition of /nl- The factors of H 2 /e in the bispectrum that are 
divided out to get A are proportional to the power spectrum which is itself oscillatory. 
This can be overcome by factoring out a number (Pk) 2 that is equal to the square of the 



value of the power spectrum that is measured when assuming scale invariance [16]. This 
gives 




«(kifc(ka)C(ks)> - i2w)'F ( > . k, ) I -f I n ^ f . (3.20) 

In the case of slow-roll and a flat spectrum this definition is equivalent to our previous one 
in Eq.( |3.18| ) (with Py. = H 2 /87r 2 e). It is then appropriate to define /nl in the same way 
as in Eq.( |3.19| ) and it is this that we discuss in the following section. Note that although 
Refs.pl, [9|, 16] define A differently, the definitions for /nl coincide. We have followed the 



convention of Ref . M for A. 



This follows directly from Eq.(l.S) if we rewrite it in k space and then evaluate the bispectrum. 
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4. Results 



4.1 Reproduction of the kink model 



To check the accuracy of our code and calculations we first set out to reproduce the 
bispectrum of the 'kink' model in Eq,( |1.3| ) |16| |. Note that cp s dictates where the kink occurs, 
c dictates the magnitude of the total shift in slope and d dictates the rate of the change in 
slope. For models where the non-gaussianity is generated as the modes cross the horizon, 
the bispectrum will be peaked over the equilateral shape (ki = k2 = k3). The authors of 
Ref.|l consider the statistic G(ki,k 2 ,k 3 )/k 1 k 2 k 3 , where / NL = _ g^.fca) . in the 

equilateral case this becomes /nl = — irl^ - m Fig-|] we show the equilateral bispectrum 
for the kink model over the full range over which it is non-zero, to be compared with Fig. 4 
of Ref.[16| which shows this over a narrower range in k. Note that although Q/k 3 ~ — /nl 
becomes as large as 10, it does not remain large over an extensive range of k. This has 
important consequences for any attempt to detect the non-gaussianity. 




Figure 1: Full equilateral bispectrum for the 'kink' model jT6|, with m = 10 6 , c = 0.0018, 
4> s = 14.81 and d = 0.022 ST) (M P = 1). The scale of k is arbitrary. 



4.2 Bispectrum of multiple inflation 

The most important, and difficult, part of the bispectrum calculation is obtaining the mode 
functions, u^. Therefore we first demonstrate that we can reproduce the power spectrum 
for multiple inflation accurately as shown in Fig.|| (to be compared with Fig. 2 of Ref . [pvfl ) . 
These spectra were calculated considering only one flat direction field, with a positive 
coupling (A > 0) to the inflaton. The parameters used in the multiple inflation potential 
flO| ) are <p = 0.01, m 2 = 0.005F 2 , A = H 2 , 7 = 1 and fi 2 = 3H 2 J35|. Note that v o 
is fixed through Eq.(p.4[) by requiring the amplitude of perturbations to be the observed 
value (~ 2.5 x 10 -9 ) after the phase transition, and is different for each value of n, the 
order of the non-renormalisable operator which lifts the flat direction potential at large 
values of the field. The time at which the phase transition occurs also varies with n. We 
emphasise that the value n = 16 is particularly favoured as it is a well known flat direction 
in the MSSM scalar potential ||]. 
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Figure 2: The power spectrum for multiple inflation |35| , reproduced for various values of the 
order n of the non-renormalisable operator in the potential (2.1). The inflation scale is chosen 



correspondingly so as to yield the same amplitude for the density perturbation on CMB scales. 



We can now confidently use our code to calculate the bispectrum of multiple inflation 
as shown in Fig.|3| Note that as k — > 10 -2 /iMpc -1 , the bispectrum begins to exhibit 
noise — this is a consequence of our not applying the regularisation scheme described in 
Ref . |Tfj|] and is only present at these scales for the reasons explained earlier. Secondly, the 
non-gaussianity, while not as large as in the kink model, is non-zero over a wider range 
of k. This can be traced to the fact that being a physical rather than a toy model, the 



features in the power spectrum generated by multiple inflation [35] are not as sharp as in 



the kink model [29]. Because the dominant term in the action is proportional to drj/dt, if 
the transition is not sharp, then the contribution to the bispectrum will not be large either. 
In multiple inflation the ip field evolves to its minimum analogous to a critically damped 
oscillator and oscillates for a period determined by the Hubble damping, which means that 
the bispectrum is non-zero for longer. 
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Figure 3: Jnl for the bispectrum of multiple inflation, with parameters that fit the WMAP 3-yr 
TT and EE power spectra f3j|. 



Due to cosmic variance (the transition occurs at large scales) and secondary non- 
gaussianities, it would be difficult to detect this signal even with a perfect estimator. We 
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have a similar concern about the bispectrum of the kink model [16] because although it has 
a larger amplitude, it is non-zero over only a decade in k making it very hard to acquire 
sufficient statistical significance for a convincing detection. 

4.2.1 The effect of varying parameters 

In Ref . |35[| , when finding the best-fit for multiple inflation to the WMAP TT and EE 
data, the only parameters varied were the amplitude of the potential, the location of the 
phase transition and the power n of the non-renormalisable term lifting the flat direction. 
All other parameters were held fixed at their most natural values, in order to maximise 



the predictive power of the model. Unlike the toy kink model [16] or other empirical 
models the parameters of multiple inflation are not free to range over any set of values 
because of its grounding in the effective field theory framework. If one simply allowed all 
parameters to vary and then performed a naive x 2 test, the physically constrained (and 
thus predictive) nature of multiple inflation would not be factored in. 

We note in passing that there is a tension between the two methods of determining 
the likelihood of a model being correct. From the fundamental physics perspective one 
can calculate the parameter values that give the best fit and then estimate how natural 
they are. From the cosmology model fitting perspective one can determine the goodness of 
fit, given the number of parameters needed to achieve it. Neither test is fully satisfactory 
because natural valued parameters could give a bad fit and a good fit could come from an 
unnatural model. A Bayesian treatment that can quantify the naturalness of a model as a 
prior probability distribution would solve this problem. The resulting Bayesian likelihood 
would incorporate both the naturalness of the underlying model and the goodness of the 
fit to the data. 

Irrespective of this issue, if multiple inflation is true, we do not know for certain what 
the values of the parameters will be. Therefore we have explored the effect of varying the 
parameters A, fi 2 and 7 in turn and show the results in Figs. [| and ||. In each plot, we 
change Vq as we vary the corresponding parameter to ensure that the power spectrum is 
always normalised to its observed value of 2.5 x 10~ 9 . For ease of comparison we also alter 
the starting point of the phase transition so that the oscillations in each figure begin at the 
same point (this does not alter the shape of the bispectrum in any way). Each plot has 
been approximately equated to spatial scales probed by the CMB observations, although 
we are not here comparing with actual data. 

The simplest parameter to consider is A which does not affect the t/j field's dynamics 
but changes the effective mass of the inflaton through -m^ = — m 2 + \ip 2 , so increasing 
A increases the magnitude of the oscillations. Therefore we expect that increasing A has 
no effect on the bispectrum except to increase its amplitude and this is just what we see 
in Fig.||. It is clear that a large bispectrum can be generated from natural values of the 
parameter A (although values above ~ 5H 2 can probably be ruled out already from the 
power spectrum). 

4.2.2 The relationship between T and Inkx 

To consider the effects of fi 2 and 7 on the bispectrum it is necessary first to derive an 
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Figure 4: Effect of parameter A on the bispectrum (n = 14, 7 = 1 and /i 2 = 3H 2 for each curve). 



identity relating the period of the oscillations over time during the phase transition to 
the period of the oscillations over k in the power spectrum and bispectrum. Eq. (|3.12| ) 
describes the time evolution of the mode functions Uf. = Vk/z of the operator £j(k, i). In 
the slow-roll approximation, z" jz = 2(aH) 2 . When this quantity is small, vj, oscillates as 



J,k(r—To) 



2k, and when this quantity is large we have Vf. = az, with a constant. 



To a good approximation, the value of the constant a is determined by equating the two 
solutions when z" / z = k 2 . This approximation will continue to hold even if z"/z 7^ 2{aH) 2 
so long as z" / z continues to grow exponentially. 

The quantity z = ay/2e can be rewritten as z = 4>'/H. Since H is almost constant we 
can write z 1 = <f>"/H and z" = <f>" /H. If there are oscillations induced on <j)' they will be 
seen more sensitively by higher derivatives of this quantity, hence z" = <j>'" /H will have a 
relative amplitude of oscillation much greater than z' = <f>'/H. This means that z" jz will 
oscillate with the same frequency as 4>' , which in turn will oscillate with the same frequency 
as if; in its minimum. 

In the slow-roll limit, a grows as e Ht and this holds even during the phase transition 
since the vacuum energy Vo remains sensibly constant. 11 Therefore, if the oscillation repeats 
itself every T units of time, then it will repeat itself every (lna)r = HT units of In a. It 
was established earlier that this means z" Jz will also repeat itself every (In o)t units of 
In a. We can thus parameterise the deviation of z" / z from 2{aH) 2 by the equation: 



z 
z 



2(aH) 2 f(lna), 



(4.1) 



where f(x) is a function with period HT. When we do the matching of the solutions inside 
and outside the horizon and equate ^- = k 2 , we get, 



— ) =k = V2aHf^, 



(4.2) 



"This will of course not hold precisely, however the relative deviation will be much smaller than the 
deviation of z" jz from its slow-roll value, due to the cj>' term in z. 
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If the period of f(x) is HT then the period of / 1//2 will be 2HT because it will repeat itself 
half as often. If we ask what the difference in k will be in the matching solution when f l l 2 
undergoes one period, from t = ti to t = i 1; with i 2 — tj. = 2T, we get, 



h = V2H 



e 2ffi2 /(2ff< 2 ) 1/2 -e 2ffil /(2ffii) 1/2 • (4.3) 



Unfortunately this result is not independent of the values of t 2 and t\ and will vary as the 
oscillations occur. However if we consider what will happen to In k during this same time 
period we get, with f n = f(2Ht n ): 

lnk 2 -lnk 1 = 2H(t 2 -t 1 ) + ^(lnf 2 -lnf 1 ). (4.4) 

Now the only time dependence comes in the difference between the function In/ at the 
two repeating points. These points were chosen because they were at the same point of 
the oscillation in / 1//2 , therefore to a good approximation In / will be equal at each point. 
This approximation will only hold exactly if the oscillations in z" / z are exactly repeating. 

It is thus expected that the mode functions will oscillate with a period over In A: 
equal to (lnfc)T = 2HT, with a small drift due to the change of the function /(In a) during 
each oscillation. This result is general to any oscillatory phase transition during inflation 
and is not specific to multiple inflation. It is a simple matter now to calculate the expected 




1 iV 10~ 3 fc(hMpc _1 ) icf z 

Figure 5: Effect of parameter fi 2 on the bispectrum (n — 16, 7 = 1 and A = H 2 for each curve). 

period, T, of the oscillations in multiple inflation. The minimum of the potential for tp 
occurs at the point, £ = (/j, 2 /nj) 1 ' ( n ~ 2 ). If we expand the potential around this point, 
with respect to the variable ifr* = — S, then the potential becomes, 

^■{n — 1) 

V(^) = V* + higher order terms. (4-5) 

It is useful to note that 7 drops out of the quadratic term and only appears in the higher 
order terms. Seen from this perspective, the equation of motion for the field ip, when in 
its minimum, is to first approximation a damped harmonic oscillator. The period, T, of a 
damped harmonic oscillator with damping term 3H and potential fi 2 (n — l) , 2 /2 is simply: 



HT = 
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From the results earlier we expect then that the period of oscillations in both the power 
spectrum and bispectrum should depend only on the ratio (fi/H) 2 and the power, n, of the 
non-renormalisable term lifting the potential of the flat direction, tp. A brief check of the 
period of P(ln/c) confirms the expectation that (ln/c)r = 2HT for the mode functions. 12 
This can be read off Fig|2] for the parameters listed in the figure. 

Therefore we expect that the result of increasing /j, 2 H 2 will be to decrease the period 
of oscillations in the bispectrum. Since this means more rapid oscillations, we also expect 
a magnification in 7/ during the phase transition, hence a magnification in the amplitude 
of the bispectrum. This is just what we see in Fig.| (th ere is also a small effect on the 
shape of the bispectrum). 




Based on this calculation we expect that changing 7 will have only a secondary effect 
on the bispectrum. Indeed in Fig.|6| the amplitude changes from ~ 5 to ~ 2 when 7 is 
changed by over an order of magnitude. By comparison, when A is similarly altered, the 
change in the bispectrum amplitude was from ~ 1 to ~ 20. The increase in the amplitude 
is because changing 7 changes the minimum of the ip field, which alters the magnitude of 
the change in the mass of the inflaton. 



4.3 The bump model of 



Finally, we point out that in Ref.[35] a second example of multiple inflation was considered 
with two flat directions which have opposite sign couplings to the inflaton and produce a 
'bump' in the primordial power spectrum. Using this model it was found that the WMAP 3- 
year data can be fitted by an Einstein de-Sitter cosmology without a cosmological constant. 
We have calculated the bispectrum for this model which has fi 2 = fi 2 = 3H 2 , 71 = 72 = 1, 
Ai = A2 = H 2 , ?7-i = 12, 77-2 = 13, and find the amplitude of the oscillations induced 



12 P(lnfc) should have period (1iiA;)t/2 because it involves the square of the mode functions Uk- To a good 
approximation, the bispectrum will have a period (lnfe)r/3; however the bispectrum is more complicated 
because the Uk are the mode functions of the £r, which are the free field operators. The bispectrum is a 
correlation function between the full operators = Ui(t,to) ' Ci(t)Ui(t,to) and there will be additional 
oscillations generated by the Ui terms. 
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by the rf term to be only ~ 0.1. Such a small value will be swamped by secondary non 
gaussianities and could never be observed. 



5. Conclusions 

There is tentative evidence that the primordial power spectrum of scalar perturbations is 
imprinted with sharp features on large scales |25], 26 1. These could have resulted from one or 



more phase transitions early on in the inflationary epoch as happens in the multiple inflation 
model ]^8| . This model is well motivated by fundamental physics (N = 1 supergravity) and 
such spectral features were predicted before WMAP provided observational indications for 
them. 



It was shown [16] that any such departure from slow-roll during inflation should also 
generate non-gaussianity. We have calculated the non-gaussianity in mutliple inflation for 
the parameter values that best fit the features in the power spectrum f3q] . This is on the 
edge of being observable but a clear detection would require a new type of bispectrum 
estimator due to its scale dependence and non factorisability jj]]. Significantly larger non- 
gaussianities can be generated during multiple inflation if the parameters are allowed to 
range over values which are technically natural (i.e. stable towards radiative corrections). 

The form of this non-gaussianity is tightly correlated with the power spectrum — the 
oscillations in the bispectrum should begin at the same multipole as the oscillations in the 
power spectrum and have two thirds of the period. There have been a number of attempts, 



to deconvolve the primordial power spectrum directly from the WMAP data [20, 21, 22 
(see also Ref.[|5l]]). It is generally found that there is a supression of power on the scale 
of the present Hubble radius, followed by a 'ringing' at medium scales. By measuring the 
period of the oscillations in the power spectrum one can predict the period (and phase) of 
oscillations in the bispectrum, in a model independent manner. 

Forthcoming measurements of CMB polarisation by Planck ought to shed light on 
whether these features in the power spectrum are systematic errors or genuine evidence 
of non-trivial dynamics during the inflationary era. The associated non-gaussian signal 
should also be detectable according to our model calculation and provide insight into the 
dynamics. 
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